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Method of filtering an image with bar-shaped structures 



The invention relates to a method of filtering an image with bar-shaped 
structures by means of Gabor filters, which are formed in the spatial domain by a two- 
dimensional Gaussian bell-shaped curve on which a cosine function is superimposed in a 
main direction. 

5 An important field of application of the invention relates to the processing of 

skin-print images, in particular fingerprint images. Starting from a gray-scale image of a 
fingerprint, as supplied by various commercially available sensors, an intensification of the 
groove structure contained is to be achieved, together with the repair of any included defects. 
In a further step, this processed image is binarized and the crucial features are then extracted 

10 from the eroded binary image. 

Gabor filters have proved their worth in the intensification and repair of (local) 
groove structures, since their principle is one of direction- and frequency-dependent 
intensification. Information concerning the main direction and main frequency is required 
locally as the input values of a local Gabor filtration, together with their local changes. The 

15 calculation of these to a sufficient quality is presupposed, and is not a constituent part of the 
present invention. 

The Gabor filtration of an image segment (e.g. a 32 x 32 tile) is normally 
undertaken with a diversion via transformation into the frequency domain (FFT = Fast 
Fourier Transformation), followed by point-by-point multiplication by the Gabor filter mask 

20 in tile size, as previously calculated from the local parameters, and then inverse 
transformation into the time and spatial domain (IFFT = Inverse Fast Fourier 
Transformation). Owing to artifacts of the transformation that occur (FFT and IFFT), the 
assembling of the individual result tiles is hereby normally insufficiently precise if this takes 
place without overlapping of adjacent tiles. Window and overlapping techniques of this kind 

25 are part of the prior art. 

In known methods, the diversion via the frequency domain is accepted owing 
to the high degree of complexity of a two-dimensional convolution in the time domain that 
would otherwise have to be accomplished. 
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It is an object of the present invention to achieve a Gabor filtration of an 
image, in particular an image of a skin print, with the least possible complexity as regards 
program implementations, wherein the filtration is to be largely locally adapted. 

This object is achieved in accordance with the invention in that the image is 
5 divided into tiles, that a predominant direction of the bar-shaped structures is determined for 
each tile and the filtration is undertaken in such a way that one tile at a time is rotated until 
the predominant direction lies at right angles to the main direction of the Gabor filter, that 
one filtration takes place in the main direction and another filtration takes place at right 
angles to this, and that the filtered tile is rotated back again. 

10 The present invention enables a local and adaptive Gabor filtration directly in 

the time domain without a diversion via the frequency domain, without a "true" two- 
dimensional convolution having to be calculated. This is enabled via the disaggregation of 
the individual two-dimensional Gabor filters into two one-dimensional filters whose vector 
product yields the Gabor filter sought, wherein the time sequence can be freely selected. 

15 However, the precise disaggregation into two one-dimensional filters is possible only if the 
wave front of the cosine oscillation and the orthogonal main axes of the superimposed 
Gaussian bell-shaped curve run in the axial direction, which is accomplished by the filter's 
fixed alignment in accordance with the invention, with the corresponding rotation of the tile. 

In contrast with the realization of the Gabor filtration in the frequency domain, 

20 Gabor filtration in the time domain is characterized by the following advantages: 

- the size of the Gabor filter can be freely selected (it does not have to 
correspond to the selected tile size), 

- the structure of the one-dimensional filters is extremely simple to realize 
with the given parameters, 

25 - the tile size is freely selectable without any extra costs arising in terms of 

time required and program memory (non-quadratic tiles are possible, not just 
dyadic, as preferably used for FFT/IFFT), 

- program is simple to realize, 

- relatively small program memory necessary, 

30 - a qualified or required accuracy (both in floating-point and fixed-point 

arithmetic) can easily be achieved, 

- no implementation of a two-dimensional FFT/IFFT for target platform is 
necessary. 
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An adaptation of the filter to the particular tile is preferably undertaken in that, 
tile by tile, for one of the filiations, a cosine oscillation with a frequency equal to the 
frequency of the structure at right angles to the predominant direction is derived, and that the 
cosine oscillation is modulated with a Gaussian bell-shaped curve, and, if applicable, tile by 

5 tile, for the other of the filtrations, the width of the Gaussian bell-shaped curve depends on 
the change in direction of the structures on the tile. 

A wider Gaussian bell-shaped curve is possible and preferable if the direction 
of the structures on the tile is subjected to only slight changes. A larger area can then be 
included in the filtration in the direction of the structures, so that irregularities that do not 

10 belong to the structures, e.g. small spots, are suppressed. In the case of contours with 
pronounced changes on the tile, i.e. strongly curved structures, however, a narrow-band 
filtration in this direction is necessary. It is further possible to set the width of the Gaussian 
bell-shaped curve in the direction of the cosine oscillation to depend on the change in 
frequency on the tile, which, however, is frequently not necessary in the case of fingerprint 

1 5 images if the tile size is matched appropriately. 

In order to reduce the implementation and calculation input further, provision 
may be made in the method in accordance with the invention for selected angles, which are 
implemented in a particular program, to be defined for the rotation, and then one of the 
defined angles that most closely accords with the rotation that is necessary per se is used for 

20 application of the filtration. Thus, for example, rotations of ±5°, ±1 0°, ±1 5°, . . . ., ±90°, which 
are minimized as regards their complexity, in particular for specific angles (±45°, ±90°), may 
be permitted. Against this background it is possible to implement the program optimally in 
terms of time for each of the selectable angles of rotation. 

Since each image element of the particular tile has to be re-calculated during 

25 the rotation in any event, the method in accordance with the invention may also be simplified 
in that, during the rotation, low-pass filtration takes place through interpolation, and/or 
binarization takes place simultaneously during the back-rotation. 

During the rotation in accordance with the invention of the tiles to be filtered, 
a greater number of image elements is inevitably included in the filtration, so that ultimately 

30 a larger tile than the resultant tile is filtered. The size of this tile depends on various 

influences, e.g. including the size of the filter. In order to filter a tile of defined size that does 
not overlap with adjacent tiles, one embodiment of the invention that is advantageous in this 
respect involves the formation of a larger tile that does overlap with the adjacent tiles and is 
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of a size of at least double the root, and, after the rotation, the larger tile is filtered in a square 
having a side length corresponding to at least double the root of the larger tile. 

A further minimization of calculation time can be made in the case of the 
method in accordance with the invention in that entries (values) lying below a threshold value 
5 and located at the edges of the one-dimensional filters are not taken into account during the 
filtration. 



The invention will be further described with reference to examples of 
10 embodiments shown in the drawings, to which, however, the invention is not restricted. 

Fig. 1 shows, schematically, individual steps of the method in accordance with 

the invention. 

Fig. 2 shows the size of a tile to be filtered at different phases of the method in 
accordance with the invention. 
1 5 Fig. 3 shows an example of a fingerprint image 

Fig. 4 shows a Gabor filter in perspective view. 



Fig. 1 shows the Gabor filtration of a tile 1 of a fingerprint image, wherein, 
20 owing to the regulations for patent drawings, the effect of the filtration itself is not shown. 
Tile 1 is one part of the fingerprint as shown in Fig. 3, and may have, for example, 32 x 32 
image elements. However, other sizes and non-quadratic, square tiles are also possible. With 
the known Gabor filtration, a two-dimensional filtration, in which the Gabor filter is aligned 
according to the bars 2 running obliquely and in a curve, would have to take place in the 
25 example shown in Fig. 1 . 

In the method in accordance with the invention, however, the mean direction 
of the bars is determined with algorithms that are known per se, whereupon tile 1 is rotated in 
such a way that the mean direction of the bars 2 is located at a right angle. The thus rotated 
tile 3 is filtered with a filter 4 comprising a cosine function modulated with a Gaussian bell- 
30 shaped curve. The frequency of the cosine function is matched, in advance, to the spatial 
frequency of the bars 2. The width of the Gaussian bell-shaped curve depends on the 
curvature of the bars and, if applicable, on the change in spatial frequency of the bars. 

In a further step, the initially one-dimensionally filtered, rotated tile is filtered 
at a right angle to the previous filtration with a Gaussian bell-shaped curve 5. The width of 
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this Gaussian bell-shaped curve depends on the mean curvature of the bars 2. Subsequently, 
the two-dimensionally filtered tile 3 ' is rotated back to its starting position. 

Fig. 2 illustrates that, in order to filter the rotated tiles, larger tiles must 
initially be derived in order to avoid artifacts at the edges. The starting point hereby is a 
5 resultant tile with 32 x 32 image elements. Owing to the rotation of the tiles, it is necessary to 
take account of a greater number of image elements. This variable is given by double the root 
as the diagonal. The filtered tile of a size 32 x 32 must be derived through back-rotation of a 
tile having a minimum size to be buffered of 50 x 50 if a bilinear interpolation is undertaken 
during the rotation. Artifacts are completely avoided if this buffer is expanded to 52 x 52 
10 image elements. 

In order to obtain a valid area of a size 52 x 52 for filtration with a filter of 1 5 
x 15 image elements, a tile to be buffered of a size 66 x 66 ultimately arises. Fig. 2 shows the 
individual sizes and rotations in undertaking the method in accordance with the invention, 
wherein it is assumed that, in order to filter a tile of 32 x 32 image elements, a tile of 46 x 46 
1 5 image elements is taken from the image. A rotation of 45° then yields a tile of 66 x 66 image 
elements, of which, however, only an area of 52 x 52 image elements can be processed with a 
filter that is 15 x 15 image elements in size. 

The broken line shows this area. From this, however, results of the filtration 
can be exploited only if all the image elements covered by the particular filter position lie 
20 within the rotated 46 x 46 image-element tile. Following back-rotation, the filtered tile of a 
size 32 x 32 emerges, which must not overlap with other result tiles of the image. 

Fig. 3 shows a black and white picture of a fingerprint image with a tile 1 , as 
processed in the example shown in Figs. 1 and 2. 

Fig. 4 shows a Gabor filter formed in one direction by a Gaussian bell-shaped 
25 curve and in the other direction, which is at right angles to it, by a cosine oscillation 
modulated by a Gaussian bell-shaped curve. 



